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(not merely the central pressure) that are considerably more subtle than might first 
be expected. 
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1. Introduction 

A useful pedagogical tool in teaching general relativity is the use of simple model- 
independent bounds on stellar structure, or more generally on compact fluid spheres. 
In the Newtonian case simple bounds are nicely summarized by Chandrasekhar, in his 
book "Stellar structure" [1], while in general relativity the most famous of the model- 
independent bounds is the Buchdahl-Bondi limit [2, 3] whereby the total mass and 
radius of any fluid body is bounded by 
2M 8 

If < 9- (1) 
Note that this is stronger than the obvious bound 2M/R < 1 that would most naively 

be expected from the absence of a black hole event horizon. 

Now the Buchdahl-Bondi limit is discussed in many of the standard textbooks (such 

as Wald "General Relativity" [4], or Weinberg "Gravitation and Cosmology" [5]) but the 

discussion can be greatly enhanced to place limits not only on the global compactness 

[= 2M/R], but also on additional features such as: 

• the "internal compactness" [= 2m(r)/r], 

• the locally measured acceleration due to gravity, 

• the metric coefficients of the spacetime, 

• various linear combinations of the pressure and density profiles, 

at all points throughout the star. The discussion below considerably extends what 
we have been able to find in the literature, but because of its elementary nature is 
particularly useful when teaching a first course in general relativity. 

Key novel results are a number of bounds on the pressure profile p(r). These 
include a relativistic extension of Chandrasekhar's bound, notable for its simplicity, and 
a number of constraints relating the pressure profile to that of the Schwarzschild interior 
geometry. We also discuss the central redshift — a quantity which is much less commonly 
discussed than the surface redshift. While the surface redshift can be measured directly 
by looking at optical spectral lines determining the central redshift would require, for 
instance, good spectral information (and theoretical modelling) regarding neutrino fluxes 
from the core. Though not a simple task, one of the payoffs would lie in thereby 
extracting information regarding the central pressure from the central redshift — an 
inequality on z c for general fluid spheres will be one of the results of this analysis. As 
always, bounds of this type £1X6 db trade-off between: 

• the ease which with the bound can be stated; 

• the ease which with the bound can be proved; 

• the ease which with the bound can be used; 

• the generality of the bound; and, 

• the strength of the bound. 

We have sought a compromise that yields useful physical information for physically 
reasonable fluid spheres. 
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2. Newtonian Stars 

For a Newtonian star the key structure equation is the equation of internal equilibrium 
[we adopt units where G = 1] 

dp p(r) m(r) m'(r) m(r) 

dr r 2 47rr 4 

Chandrasekhar [1] , following on and adapting earlier classical work, uses this to establish 
several theorems of which one stands out as the most important for the present purposes: 

Theorem 1 Provided the average density p{r) = m(r)/[(47r/3)r 3 ] is nonincreasing as 
one moves outward in the star, then so is the quantity 

3 m(r) 2 



p(r) + 



8tt 



Proof: 



dr 



p(r) + 



3 m(r) 2 
8^ r 4 



dp 

dr ' 47T 

1 m(r) 

~ 2^ r 4 

1 m(r) d 

2n r dr 

2m(r) d 

= 3~ 
< 0. 



3 m(r) m'{r) 3 m(r) s 



2vr 



m'{r-) 



3m(r) 



m r 



dr 



(3) 

(4) 

(5) 

(6) 

(7) 
(8) 



Corollary 1 In a Newtonian star where the average density is nonincreasing outwards 

(9) 



, , 3 m(r) 2 3 M 2 
Pc > p(r) + — > 



8tt 



r-4 - 



8tt i? 4 



The net result of this theorem is to give robust and useful bounds on the pressure and 
density distribution inside a Newtonian star with only minimal input information. We 
will seek similar robust and simple results for general relativistic stars. 



3. TOV Stars 

The key structure equation for a relativistic star is the Tolman-Oppenheimer-Volkoff 
[TOV] equation 

Mr) = [p(r) +p(r)} [m(r) + 47rp(r) r 3 ] 
dr r 2 [l -2m(r)/r] ' 1 ] 

The most obvious point to make is that because density, pressure, and mass are all 
assumed positive in any reasonable star we now have the inequality 



dp(r) p(r) m(r) m'{r) m{r) 



dr 



47rr 4 



Consequently the Newtonian theorem and its corollary still holds in a relativistic star 
- the only modification is in the proof of the theorem; where the equality on going 
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from line 1 to line 2 is now an inequality. Of course this is a very weak-constraint on the 
structure of a relativistic star, and much better constraints are available. A particularly 
simple constraint is based on the fact that for a relativistic star the positivity of pressure 
implies 

dp(r) ^ p{r) m(r) m'(r) m(r) 

dr ~~ r 2 [l — 2m{r)/r] 47rr 4 [l — 2m{r)/r] 
Note that by dropping both linear and quadratic pressure terms on the RHS we are 
implicitly ignoring the "regeneration of pressure" [6]; consequently the bound we derive 
will be relatively weak but has the advantage of an elementary proof. 

Theorem 2 Provided the average density p(r) = m(r)/[(47r/3)r 3 ] is nonincreasing as 
one moves outward in a relativistic star, then so is the quantity 



p(r) - 



3 m(r) 



16n r 3 



In 



1 - 



2m{r) 



Proof: 



3 m(r) 



16tt 
dp 
dr 



In 



2m(r) 



(13) 



(14) 



3 d 
167T dr 
3 

+ 



m(r) 
m(r) 



In 



2m(r) 



r 
d 



< 



1 m(r) 



16n : 
1 



2m(r) 



3 1 — 2m(r)/r dr 

3m(r) 



8n r 4 1 — 2m(r)/r 
1 m(r) 1 



m r — 



8tt 
< 0. 



2m{r)/r dr 



(15) 

(16) 

(17) 
(18) 



Corollary 2 In a relativistic star, provided the average density is nonincreasing as one 
moves outward 



Pc > p{r) - 



3 m(r) 
16-7T r 3 



In 



1 - 



2m(r) 



> 



3 M 
16tt R 1 



In 



1 - 



2M 



(19) 



Comment: 

Note that as 2M/R — > one recovers the Newtonian result, and that as 2M/R — > 1 the 
central pressure is forced to infinity. This argument is not good enough to reproduce 
the Buchdahl-Bondi limit, but it is elementary, and it is sufficiently good to show 
the runaway explosion of central pressure for sufficiently compact objects. Moreover, 
because in deriving this bound we have not used the "regeneration of pressure" effect [6], 
it follows that the occurrence of unbounded central pressures in general relativity is not 
intrinsically a regeneration effect. 

To now obtain the Buchdahl-Bondi bound, and its extensions, the best strategy 
is to develop a number of comparison theorems that provide inequalities relating the 
generic relativistic fluid sphere and the interior Schwarzschild solution. 
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Let us adopt coordinates to write the general interior geometry in the form 

„2 _ /-/ \2 1.2 



dr 2 



oV = -C(r) 2 dt 2 + . . . + r 2 [d9 2 + sin 2 9 d0 2 ] . (20) 

1 - 2m(r)/r 1 J 

Equivalent ly 

ds 2 = - exp -2 / g(f) df dt 2 + — — + r 2 \d9 2 + sin 2 9 d0 2 l . (21) 

J r 1 — 2m(r)/r L J 

With this choice g(r) is positive for an inward gravitational attraction and equals the 
locally measured acceleration due to gravity. The Einstein equations yield 

Svrp = G u = 2m'(r)/r 2 ; (22) 

8n P = G ff = 2 j ^ (r) [1 - 2m(r)/r] - ^} ; (23) 

m'(r)[l + r^(r)] m(r)[l — rg(r) — 2r 2 g'(r) — 2r 2 g(r) 2 } 
8irp = G §§ = + 

| g(r)[l + rg(r)} ^ 

The first of these integrates to 

//Mr) = / Anp(f) f 2 df, (25) 
Jo 

which justifies the choices of notation g rr = 1/[1 — 2m{r)/r\. The second equation can 
be algebraically rearranged to yield 

= m(x) + 4?rp(r)r 3 
9{T) r 2 [1 -2m{r)/r\ 1 j 

The third equation is quite messy — however in view of the Bianchi identities it may be 
"traded off" for the covariant conservation of stress energy which [for isotropic pressures] 
takes the rather simple form 

p\r) = ~[p(r) +p{r)\ g{r) (27) 

Eliminating g{r) leads to the TOV equation as previously given. 

If we choose to work in terms of £(r) instead of g(r), and we shall soon see that 
this is sometimes useful, the only significant change in appearance [not substance] is in 
the equation 

_ C(r)"[l - 2m{r)/r\ ('(r) [1 - m(r)/r - m'{r)\ [m(r)/r]' 

~ (F) + rW) (28) 

For the G?? equation there is very little change 

„ n ( C'(r) [1 — 2m(r)/r] m(r)] . . 



Bounds on the interior geometry and pressure profile of static fluid spheres 



6 



5. Interior Schwarzschild 



The interior Schwarzschild solution is a specific solution of the Einstein equations 
corresponding to a star whose density is constant throughout space. (Despite the 
common misperception, this does not mean the star is "incompressible", see Misner, 
Thorne, and Wheeler for details [7], see esp pp. 609 ff.) If we write the metric as 

dr 2 



dsi = -UrY dt' + 



+ r 2 [d6 2 + sin 2 9 d(j) 2 ] , 



(30) 



1 — 2m»(r)/r 

then solving the Einstein equations for p = p* yields the standard textbook results [4, 
5, 7]: 



m*(r) 



— p*r 3 = M —, 



C*(r) = - [3^1 - 2M/R - y/1 - 2m*(r)/r 



and 



P*{r) = 



y/l-2m*(r)/r-^l-2M/R 



(31) 
(32) 

(33) 



3y/l -2M/R- y/l-2m*{r)/r' 

Here and henceforth a subscript or superscript star on a quantity refers to the 
Schwarzschild interior solution. The central pressure is 

1 - Jl - 2M/R . . 

p* = p* , V 34 
3^1 -2M/R- 1 

which certainly diverges as 2M/R ^8/9; but this is not yet enough to establish the 
Buchdahl-Bondi bound — to do so we need to develop theorems relating the general 
interior geometry (for an object with the same mass M and radius R) to the specific 
Schwarzschild interior geometry. At the center of the interior Schwarzschild geometry 
we also have 

C = \ [3 Vl - 2M/i? - l] (35) 
which is related to the central redshift by 



a" 1 



p* 



As useful algebraic results note that 

P* 



p*(r)C*(r) = y \yi-2m*{r)/r-y/l-2M/R 



and 



[P*+P*(r)} C*(r) = p^l-2M/R, 
[p* + 3p*(r)] C*( r ) = P*a/1 - 2m*(r)/r. 



(36) 

(37) 
(38) 

(39) 
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6. Comparison theorems for the spacetime geometry 

Now consider a generic self-gravitating relativistic fluid sphere 
Lemma 1 // the average density p is nonincreasing outwards, then 

m(r) > m*(r) = M — . (40) 

Proof: 

By definition and hypothesis 

, , 3m(r) 3M 

Lemma 2 T/ie average density p is nonincreasing outwards if and only if for all r 

p(r) > p(r). (42) 

Proof: 

By definition and hypothesis 



dp(r) d 
dr dr 



3m(r) 
4nr 3 



I W) - p(r)} . (43) 



Note that if the unaveraged density p(r) is nonincreasing outwards then certainly p > p 
and so the average density p(r) is also nonincreasing outwards; the converse however does 
not hold. A nonincreasing average density does not necessarily imply a nonincreasing 
density and is a genuinely weaker constraint. From the first lemma it trivially follows 
that: 

Lemma 3 If the average density p is nonincreasing outwards, then 

9rr(r) > g* rr (r). (44) 

and 

m{r)/r > m*(r)/r. (45) 
Much more subtle is the following result regarding the tt metric component: 
Theorem 3 If the average density p is nonincreasing outwards, then 

C(r) < C.(r). (46) 

Equivalently 

\9u(r)\ < \g* tt (r)\. (47) 

Proof: 

Obtaining this result is somewhat tedious but relatively straightforward ■ an 
adaptation of the discussion in Weinberg [5] or Wald [4] suffices. With a little more 
work, you could also adapt the discussion in the original 1959 Buchdahl paper [2]. 
Consider 

= [G,r - G §§ ] r 3 C(r) (48) 
= — r 3 ("(r)[l — 2m(r)/r] + C'( r )[ r2 + r 2 m'(r) — 3rm(r)] + ((r)[rm'(r) — 3m(r)]. 
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Rearranging 



dr 



r y r 



C(r) 



a/1 - 2m(r)/r dr 



m(r) 



So, as long as the average density is nonincreasing as we go outward 

< 0. 



d 
dr 



r 



2m(r) 



(49) 



(50) 



At the surface of the star £(r) must match smoothly onto the exterior Schwarzschild 
geometry, therefore 

M 



C(R) 



2M 



1 



Combining (50) with this boundary condition implies [Vr G (0, i?)] 

m*(r) 1 



CV) > 



M r 



1 



R 3 a/1 - 2m{r)/r r 2 ' ^/l - 2m(r)/r 



Integrate from r to -R, [no£ from to R], we have 

M rR 
C(R) - C(r) > ^ 



a/1 - 2m(f)/r 



dr. 



But in the RHS, m(r) > m*(r), so 

M ^ 
C(i?) " C(0 > ^ 



a/1- 2m, (f)/ 



dr 



r 



""^X 7i 



2Mr- 2 jR? 



dr. 



This integral is now do-able in closed form, with the result 

R 



C(R) - C(0 > 



1 



1 



2Mf 2 



Rearranging, and applying the boundary condition at R, we have 

1 



C(0< 



3W1 



2M 



-A 1 



2Mr 2 



(51) 



(52) 



(53) 



(54) 



(55) 



(56) 



The RHS is now recognizable as C*( r )) the quantity appropriate to the Schwarzschild 

interior solution, so the theorem is proved. 

Comment: 

- Note what has happened here: Both of the [physically nontrivial] metric components, 
gu and g rr , have been bounded in terms of what their values would have been for a 
Schwarzschild interior solution of the same mass and radius. This now is more than 
enough to deduce the Buchdahl-Bondi bound. Since ((r) < C*( r )j and since C*( r ) ~~ *■ 
for some finite positive r whenever 2M/R > 8/9, we deduce that any star [in which the 
average density is nonincreasing outwards] must likewise satisfy the same bound. More 
formally: 
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Corollary 3 If the average density p is nonincreasing outwards, then 
2M 8 



R ~ 9 
Proof: 

If the center of the star is to be regular then we must have 



(57) 



< Ce < C 



(58) 



Comment: 

- It might be tempting to conclude that for all r 



III 



2m{r) ^ 8 
r ~ 9 



??? 



(59) 



Although true, this conclusion cannot be drawn from the arguments so far presented. 

- Note also what has not happened here. We cannot [at this stage] deduce any bound 
on the pressures. Indeed it might be tempting to assert 



III P(r) > p*(r) ??? 
but such an assertion would actually be false. To see this compare 



and 



dp 
dr 

dp* 
dr 



p{R) 



M 



1 



R 



R 2 1-2M/R' 



P* 



M 



R 2 1 



> 



dp 
dr 



(60) 



(61) 



(62) 



R 



, 2M/R 

Thus sufficiently near the surface we must actually have p(r) < p*(r); though at the 
center (as we shall soon see) p(0) > p*(0). Thus general bounds on the central pressure 
and pressure profile are trickier to establish (which is why, regarding the central pressure, 
Weinberg resorts to "it can be shown that", and Wald remains silent). 

As well as the metric components, there is a similar constraint on the locally 
measured acceleration due to gravity detected by static observers. 



Lemma 4 // the average density p is nonincreasing outwards, then 
g(r) > g*(r). 

Proof: 

To see this use the bound (52) to obtain 



(63) 



C(r) > 



m*(r) 



1 



> 



m*(r) 



1 



' 2 a/1 - 2m(r)/r r 2 ^/l - 2m*(r)/r 
and then combine this, the definition of g, and (46), to obtain 



g{ r )= m>m 



g*{r). 



(64) 
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2m(r) 8 
~r ~ 9 



y/l + Q7rp(r)r 2 1 + Q7ip(r)r 2 



(65) 



Proof: 

This particular result is derived, for instance, in Wald, "General Relativity" [4]. To 
obtain a bound on 2m(r)/r modify the preceding argument by now considering the 
region f G (0, r) with r < R. Then, following the discussion in Wald or using equation 
(50) above we have for r < r 



The RHS evaluates to 

1 [ 2m(r) . . 1 _. . m(r) + 4irp(r)r 3 Anpir) + m(r)/r 3 > . . 

rV r ryi-2m(r)/r yjl-2m{r)/r 
Thus 



C'(f) > C(r) 



47rp(r) + m(r)/r 3 



a/1 — 2m(r)/r a/1 — 2m(r)/r 
It is now sufficient to integrate f from zero to r: 



That is 



CW-C(0)>C(r 



C(0) < C(0 1 - 



^47rp(r) + m(r)/r 3 

yT — 2m(r)/r Jo a/1 — 2m(f)/f 

4np(r) + m(r)/r 3 



dr. 



dr. 



(66) 



(67) 



(68) 



(69) 



(70) 



a/1 — 2m(r)/r Jo a/1 — 2m(r)/r 

But because £(0) must be positive in a star whose geometry remains regular all the way 
to the center, the quantity in braces must be positive, and [independent of the unknown 
value of C( r )] we have the inequality 



I 2m (r) ^ 



, . m(r) 
47rp(r) + 



o a/1 — 2m(r)/r 



dr. 



(71) 



But, because the average density is nonincreasing outwards, our first lemma implies 

^3 



m(r) > m(r) — . 



This implies 



2m(r) 



> 



A-np{r-) + 



m r 



o a/1 — 2m(r)r 2 /r 3 



dr. 



The integral is now elementary, with the result 

m(r 



2m (r) ^ 



47rp(r) + 



2m(r) 



2m(r) 



(72) 



(73) 



(74) 
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Consequently 

^/l-^> Anp(ry + 

Equivalent ly 



2 , m ( r ) 



(75) 



1 - ^ + Jl - ^ > 47 rp(r)r 2 + (76) 
r y r • r 

It is now a matter of tedious but straightforward algebra to rearrange this inequality 
into the desired form. Finally, because for positive x one has a/1 + x < 1 + x we deduce: 

Corollary 4 // £/ie average density p is nonincreasing outwards, then 

r ~ r ~ 9 ' 1 J 

Proof: 

Simply use the fact that p(r) is positive inside the star, and the previously derived lower 
bound on the compactness. If you prefer you can write this as 

2M 4 <^< 5. (78) 

Comment: 

Since 2m(r)/r is zero at the center, and equals 2M/R at the surface it is tempting to 
suppose that 2m(r)/r might be monotone increasing. In general it is not, and many 
equations of state are known for which 2m(r)/r develops damped oscillatory behaviour 
as a function of r [6, 8]. 

Corollary 5 If the average density p is nonincreasing outwards, then 



PM<^-^ + \A-^. (79) 
Comment: 

This bound is somewhat unusual in that it provides an upper bound on the pressure. 
It is less useful than one might imagine since it gives no information about the central 
pressure. 

8. Comparison theorems for the pressure profile 

General and stringent theorems regarding the central pressure and pressure profile 
[rather than metric coefficients and compactness] are relatively more cumbersome. 
With the tools we have at hand, an easy result is this: 

Theorem 5 // the average density p is nonincreasing outwards, then 
[p(r)+3p(r)]C(r) > ^ = [p« + 3 p«(r)]C«(r) ^ 
a/1 — 2m(r)/r a/1 — 2m»(r)/r 
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Proof: 

First, by the geometric comparison theorems of the previous section, in particular 
equation (52), 

(r) = C'jr) > Mr 1 = Anp.r 1 

9{T> ((r) ~ R 3 y/l-2m(r)/r{(r) 3 y/1 -2m(r)/r £(r)' 
Second, by the ff Einstein equation 

1 j g(r)[l — 2m(r)/r] ml 1 g(r)[l — 2m(r)/r] p(r) 
P= A^{ r 4^ r 3~" ( ' 

Combine. We have 



C(r) 

For the last equality use the properties of the interior Schwarzschild solution, as 
embodied in equation (39). 

Corollary 6 If the average density p is nonincreasing outwards, then 

pc + 3 p c 
For the central redshift 

z c < Pc + 3 Pc-P* . (85 ) 
p* 

Corollary 7 If the average density p is nonincreasing outwards, then 
p{r) + 3 p(r) > p. = p* +3p,(r) 
y/\ - 2m{r)/r ~ C*(0 ^1 - 2m,(r)/r' 

Proof: 

Apply the inequality £(r) < C*( r ) to the preceding theorem. 

Corollary 8 If the average density p is nonincreasing outwards, then 

Pc + 3p c > p* + 3p* c . (87) 

Comment: 

While somewhat crude, this bound has the benefit of being both elementary and 
guaranteeing that the central value of p + 3p diverges at or before one reaches the 
Buchdahl-Bondi limit. Note that this does not permit us to place a bound on p c itself, 
only on the combination p c + 3p c . 

The first bound we place on the pressure profile itself is this: 

Theorem 6 // the average density p is nonincreasing outwards, then 

p*2y/\- 2m{r)/r 2Mr 2 /R* - 3^1 - 2M/R 

p(r) > — = — - — . (88) 

3 3^1 - 2M/R - y/1 - 2Mr 2 /R 3 
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Proof: 

Start from the equation p' = — [p + p]g = — [p + p\C/( and re- write it in the form 

[K]' = -p£. 

Integrate from r to R, then 



P(r) ({r) 



p(r) ('(f) dr. 



But because of our earlier inequality on C'( r )) embodied in equation (52), 



p(r) C(r) > / p[r 



M 
R* 

P*_ 
3 



R 3 y/l - 2m(f)/f 
m'(f) 



dr 



R 



Anr^l — 2m(r)/r 
m(f) 



dr 



m r 



1 



= y \ - a/1 - 2m(r)/r 



R 



+ 



a/1 — 2m(r)/r 
R m(f)/f 2 



dr 



dr 



y/1 - 2m(r)/r 

Again using the fact that the average density is nonincreasing outwards 



p(r) C(r) > <j - a/1 - 2m(r)/r 



R 



+ 



R 



Mf/R? 



dr 



(89) 
(90) 

(91) 
(92) 
(93) 
(94) 

(95) 



/ l-2Mf 2 /R 3 

The remaining integral is again elementary, with the consequence 
P(r) C(r) > j | -2 a/1 - 2M/R + 2a/1 - 2m(r)/r - yf\ - 2M/R + a/1 - 2Mr 2 /i? 3 } 

(96) 



= ^ {2a/ 1 - 2m(r)/r + yf\ - 2Mr 2 /R 3 - 3y/l - 2M/r} . 



(97) 



Finally, since we have already established ((r) < C*( r ), inserting the explicit form of 
C*(r) completes the proof. 

Corollary 9 // the average density p is nonincreasing outwards, then 



((r)p(r) - ^a/1 - 2m(r)/r > C*(^)p*(^) - tta/1 - 2m 3f (r)/r (98) 
Proof: 

This is a slight strengthening of the previous theorem obtained by not invoking 
C( r ) ^ C*( r ) & t the last step. In this version of the result everything on the LHS 
refers to the body we are investigating and everything on the RHS refers to the interior 
Schwarzschild geometry we are using to compare it to. 



p* 



Corollary 10 // the average density p is nonincreasing outwards, then 



p(r) > p*{r) + y 



y/j - 2m(r)/r - a/1 - 2m t (r)/r 



(99) 
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Proof: 

This is simply a re-writing of the previous theorem to make it clear that the generic 
pressure bound consists of the interior Schwarzschild result plus a "correction" term. 
The correction term vanishes at both R and at 0, and is negative between these locations. 
Note that the bound is sharp in the sense that the interior Schwarzschild solution 
saturates the bound. 
Comment: 

Note that to bound the pressure profile p(r) you need information regarding the local 
compactness m(r)/r, which is difficult to measure observationally. This contrasts with 
most of the other bounds derived in this article which depend on parameters that can 
be measured from the stellar exterior. 
Comment: 

If we take this pressure bound and consider it near the surface of the star then we find 
p(R) > 0, (100) 

and 

dp 
dr 

But, because of the boundary conditions at the surface of the star, these bounds are 
actually saturated by any fluid sphere, and cannot be improved. 



> -p(R) -777 — — — --. (101) 



Corollary 11 If the average density p is nonincreasing outwards, then 

Pc > Pi- (102) 

Comment: 

Thus the interior Schwarzschild solution really does provide a lower bound on the central 
density even if it cannot be used [in "uncorrected" form] elsewhere within the star. 
Consequently, in the Buchdahl-Bondi limit 2M/R ^8/9 not only does ( c — > 0, but we 
also have p c — > oo [as intuitively expected]. 

We now consider several additional bounds that are derived in slightly different 
fashion. 

Theorem 7 If the density p [not the average density] is nonincreasing outwards, then 
in terms of the surface density p s : 

Jl - 2M/R Jl - 2M/R , . 

p(r) + P {r) > Ps V c(r) > Ps V Ur) ■ (103) 

Proof: 

Start from the equation p' — — [p + p] g = — [p + p\('/( and re- write it in the form 

[(p + p)C]' = VC<o (104) 

Integrate from r to R, then 

p s y/l - 2M/R - \p{r) + p(r)\ C(r) < 0. (105) 
Rearrange, and use the previous inequality for C( r )- 
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Comment: 

This is a relatively weak bound because it involves the surface density p s . On the other 
hand, the RHS of the bound is now expressed in terms of quantities that are in principle 
accessible via external observation. 



Corollary 12 If the density p is nonincreasing outwards, then 

Jl - 2M/R 2J\ - 2M/R 
Pc + pc> ps ji = Ps 



C rS 3 V /1-2M/R-1 

The theorem above can be considerably strengthened by noting that p' — — [p + p]g 
— [p +p]C/C implies the exact general result 



(106) 



C(r) = C c exp 



Jo 



dp/dr 
p(f) + p(f) 



dr 



(107) 



But, assuming the density is nonincreasing outwards, we have p c > p(r) > p(r), whence 



C(r) > C c exp 

Similarly 

C(^) < C c exp 
We summarize this as a theorem: 



dp/dr 
o Pc + p{r) 

dp/dr 



dr 



Cc 



o p(r)+p(f) 



dr 



c 



Pc+Pc 

Pc + p{r) ' 
p(r) +p c 



p(r) +p(r) 



(108) 



(109) 



Theorem 8 // the density p [not the average density] is nonincreasing outwards, then 



C(r) > Cc 



Pc+Pc 



and 



C(0 < Cc 



p(r) + Pc 
p(r)+p(r)' 



(110) 



pc+pij-y 

If we now evaluate this at the surface of the body: 

Corollary 13 // the density p [not the average density] is nonincreasing outwards, then 
y/1 - 2M/R > C c ^±^; and y/l - 2M/R < ( c ^i±^£. (m) 

Pc ' Ps 

A closely related result is obtained by similar bounds applied to the exact general result 



C(r) = v 7 ! - 2M/R exp 



R 



dp/dr 



dr 



(112) 



p(f)+p(f) 

Without repeating the details we simply quote the results. 

Theorem 9 // the density p [not the average density] is nonincreasing outwards, then 



and 



C(r) > y/1 - 2M/R 



C(r) < v 7 ! - 2M/R 



p s +p(r)' 



p(r) 



(113) 



(114) 



p(r) +p(r)' 

These constraints bound the metric components in terms of the pressure and density 
profiles, and vice versa. A qualitatively different type of constraint is provided by the 
following result. 
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Theorem 10 If the sum of the pressure and density is positive 



y/1 - 2m(r)/r > C(r) > ( c y/l - 2m(r)/r 
Proof: Consider 



dr 



C(r) 



y/l — 2m(r)/r 
C(r) 

[1 - 2m(r)/r] 3 / 2 
C(r) 



C(r) 



a/1 — 2m(r)/r 
m + 4iipr 3 4npr 



9 + 



[m{r)/r\' 
1 — 2m{r)/r 



m 



Therefore 



47rr 



[1 - 2m{r)/rf/ 2 
C(0 



[p + p] >o. 



< 1. 



16 
(115) 

(116) 

(117) 
(118) 

(119) 



y/l - 2m(r)/r 
Comment: 

Note that this bound makes no reference to the interior Schwarzschild solution for 

comparison purposes. 

Comment: 

The condition p + p > is the so-called null energy condition [NEC] - - it is certainly 
satisfied in bulk nuclear matter, though there are exotic situations (typically quantum) 
in which the NEC is violated [9]. For the purposes of the present article the NEC is 
automatically assumed in view of our basic hypotheses that pressure and density are 
both individually positive within the fluid body. 



Corollary 14 If the sum of the pressure and density is positive 
£<\gtt{r)\g rr {r)<l. 



(120) 



We can now tighten this result by keeping track of the terms we know to be positive a 
little bit longer. Re- write equation (118) in the form 



Whence 



y/l - 2m{r)/r d 
C(r) dr 

C(0 



C(0 



y/1 - 2m(r)/r 



4irr 



p + p 



exp 



-r 

J r 



1 — 2m{r)/r 
p(r) + p(f) 



Aur ' v ' ' l . \ ' df 
1 — 2m{r)/r 



y/1 - 2m{r)/r 

To re-write the RHS, we use the TOV equation written in the form 
dp Anr(p + p)(p + 3p) 
dr 3(1 — 2m{r)/r) 

to obtain 

C(r) [ f R 3 dp/dr 



y/l - 2m{r)/r 



exp 



p(r) + 3p(f) 



dr 



(121) 



(122) 



(123) 



(124) 
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But now if we use our standard assumption that the average density is nonincreasing 
outwards we have first, p > p*. Therefore l/(p + 3p) < l/(p* + 3p), and because dp/dr 
is negative, p'/ (p + 3p) > p'/ (p* + 3p). Then 

C(r) \ f R 3 dp/dr 



a/1 — 2m(r)/r 



> 



exp 



dr 



p* 

p* + 3p(r) 



(125) 



p* + 3p(r) 

Conversely, since f e {r,R), we have p(f) < p(r). Therefore l/[p(f) + 3p(f)] > 
l/[p(r) + 3p(f)], and because dp/dr is negative, p'/[p(r) + 3p(f)] < p'/[p(r) + 3p(f)]. 
Then performing the integration yields 



CM 



< 



exp 



R 



3 dp/dr 
p(r) + 3p(f" 



dr 



p(r) 



p(r) + 3p(r) ' 



a/1 — 2m(r)/r 
and we have our penultimate theorem: 

Theorem 11 If the average density p is nonincreasing outwards, then 



(126) 



C(r)[p* + 3p(r)] 



> p*; 



and 



C(r)[p(r) + 3p(r)] 



< p(r). 



(127) 



y/l - 2m{r)/r ' y/1 - 2m{r)/r 

Note that the RHS of the first inequality is independent of location within the gravitating 
body. Note that this result is definitely more restrictive than equation (80), which 
was obtained by somewhat simpler techniques. This theorem has several interesting 
corollaries, easily derived using the comparison bounds on £(r) and m(r). It is easiest 
to first write 



Corollary 15 If the average density p is nonincreasing outwards, then 



C(r) [p* + 3p(r)] 



> P* 



C*( r ) [p* + 3p*(0] 



(128) 



y/1 - 2m(r)/r y/l-2m*(r)/r 
Then, since ((r) < (*(r): 

Corollary 16 // the average density p is nonincreasing outwards, then 

P* + M r ) > P* + 3p*(r) , 12Q . 
y/l - 2m(r)/r ~ y/l - 2m»(r)/r ' 

Similarly, using the fact that m(r) > m^r): 

Corollary 17 If the average density p is nonincreasing outwards, then 

- P(r) > - P * {r) . (130) 
\Jl — 2m{r)jr y 1 — 2m»(r)/r 

Though the derivation was somewhat complex and lengthy, this last result is by far the 
cleanest looking bound we have derived on the pressure profile. 
Finally, for the central redshift we have: 

Corollary 18 If the average density p is nonincreasing outwards, then 

P* J > ^ Pc 



Cc> 



P* + 3p c ' 



and 



Pc + 3p c ' 



(131) 
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so that 



and 



3pc , , 3p c 

— < z c < — , (132) 

Pc P* 



P* z c , , Pc /ioo\ 

< Pc < — =— • (133) 



3 - JU - 3 

Finally, for our last theorem we start from the exact result 

3 dp/dr 



- C(r) =C c exp 
- 2m(r)/r 



o p(*9 + 3p(r) 



dr 



(134) 



and again bound p c > p(r) > p(r). Without repeating the details is is now easy to see 
that: 

Theorem 12 // the average density p is nonincreasing outwards, then 

CW >Cc^±^, (135) 



v/l - 2m(r)/r ~ Pc + 3p(r) ' 

and 



c(r) < c f: )+3 ^v (ise) 



v/l - 2m(r)/r ~ p(r) + 3p(r) ' 
9. Conclusion 

We have extended classic work by Buchdahl [2] and Bondi [3] to derive several model- 
independent limits on stellar structure that go beyond what is commonly found in 
the literature. (Several important papers are those of Kovetz [10], Islam [11], and 
Forrester [12]; despite the wide variety of inequalities derived in those papers the 
constraints on the pressure profile derived above seem to be new.) Many of the theorems 
in this present article are of the form "If the average density is nonincreasing outwards, 
then generic stellar structure is in many ways bounded by the interior Schwarzschild 
solution of the same mass and radius". While this is not unexpected, classic results 
typically refer to global aspects such as the overall compactness [= 2M/R] or the 
central pressure p c — in this article we have seen how to extend this to local information 
regarding the metric coefficients and, local compactness, and pressure profile at all points 
interior to the star. The bounds we derive are sharp in the sense that they are saturated 
by the interior Schwarzschild geometry, so that in this sense they are the "best possible" . 

Interest in these issues arose while teaching a course in general relativity and noting 
that while many textbooks discuss the Buchdahl-Bondi bound, most did not do so by 
name, and no extant textbook [or, so far as we are aware, the technical literature] covers 
all the inequalities derived herein. 
Further issues that might bear looking at: 

— Are there any comparable bounds under more stringent hypotheses? 

(For Newtonian stars making the stronger assumption that p itself is nonincreasing 

yields somewhat better bounds on the pressure profile [1]. We have used this hypothesis 
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for some of our theorems but suspect that even better bounds may be possible.) 

— Are there any comparable nice upper bounds on the pressure? 

(We have already seen one upper bound on the pressure profile which unfortunately 
never gives any information about the central pressure. For some progress on this 
question see Kovetz [10].) 

— Are there similar bounds for solid objects such as planets? 

(This leads into the subject of "anisotropic stars" for which some of the best bounds 
available seem to be those of Guven and O'Murchadha. See [13] and references therein.) 

— Are there similar bounds for a liquid star enclosed in a solid crust? 

(This would require some hybrid techniques — a fluid interior with an anisotropic crust.) 

In closing, we reiterate that while the techniques used to derive these bounds are 
"elementary" in the sense that they can be explained in a introductory course on general 
relativity, the results are rather non-trivial. 
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